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Abstract 

Schrodinger bound-state problem in D dimensions is considered for a set of central 
polynomial potentials containing 2q arbitrary coupling constants. Its polynomial 
(harmonic-oscillator-like, quasi-exact, terminating) bound-state solutions of degree 
N are sought at an (g+1)— plet of exceptional couplings/energies, the values of which 
comply with (the same number of) termination conditions. We revealed certain hid- 
den regularities in these coupled polynomial equations and in their roots. A par- 
ticularly impressive simplification of their pattern occurred at the very large spatial 
dimensions D ^ 1 where all the "multi-spectra" of exceptional couplings/energies 
proved equidistant. In this way, one generalizes one of the key features of the ele- 
mentary harmonic oscillators to (presumably, all) non- vanishing integers g > 0. 
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1 Introduction: quasi-exact terminating solutions 

The never-ending story of the search for exact bound-state solutions started with 
the very emergence of quantum mechanics. Its part which pays attention to the 
polynomial central potentials V{r) in the ordinary differential "radial" Schrodinger 
equation 

-^"ir) + ^^^^jjj{r) + V(r)jlj(r)^E^Ij(r), ^£^^2(0,00) (1) 

is not much younger. Indeed, the elementary nineteen-century mathematics proves 
sufficient for the construction of ■?/'(r) from (1) in analytic form with, say, a power- 
series ansatz for components A{r) and B{r) in 

^(r) = r^+i A{r) e^^"-^ . (2) 

The most common harmonic-oscillator model V''^^\r) = u;^r^ provides a particu- 
larly appealing illustration of such an approach because the semiclassical exponent 
itcT^ in eq. (2) describes the correct asymptotic decrease of ^''^'^'^^(r) 
while the Taylor series for A^^'^^ (r) degenerates to a polynomial as well. 

A broad family of polynomial potentials admits a similar specification of their 
"asymptotically optimal" polynomial exponents B{r). Vice versa, for all the "canon- 
ical" polynomial WKB-like exponents 

S(^^^)(r) = iaor^ + ^a.r^ + ... + ^".-^^^^ (3) 

and for all the "canonical" power-series choices of the ansatz (2), 

00 

V'(r) = E exp i-BwKBir)] (4) 

n=0 

potentials may be polynomials with 2g -|- 1 arbitrary couplings, 

V{r) = {r)^gor^ + gir^ + ... + g2g r^''^^ = [O^"^ (r)]^^ + (r) . (5) 

The first, asymptotically dominating auxiliary factor 

n^i) (r) = ao + air^ + . . . + r^'' (6) 

is determined precisely by the plet of the WKB-related free parameters while 

5(9) (r) = Gor" + dr^ + . . . + G^-ir^^ 

carries just the asymptotically less relevant information about the full force V{r) at 
any q > 1. All the relevant details may be found in our older review of the related, 
so called Hill-determinant bound-state method [1] . 



1 



Due to the one-to-one correspondence g2q = a^^, g2q-i = g2q-i{<^q, <^q-i) = 
2 ctq-i Qiq, . . . etc (or, in opposite direction, aq = ^Jg^q > 0, aq-i = g2q-i/{2aq) 
etc.), we may work with both the old and new couphngs. Moreover, using the trivial 
changes of variables in our differential equation (1) (r'^ = x etc, with all details de- 
scribed again thoroughly in the above-mentioned review [1]), the canonical potential 
(5) generates the whole series of its mathematical equivalents, 

(x) ^fox-^ + hx + f2X^ + ... + /2, /2, > , (7) 

M (Z) = ho Z-^/^ + hiZ-' + ... + h2q-l Z'^-^'^ + h2q z'^-^ k^q >0 , (8) 

etc. Thus, the well known one-to-one mapping between harmonic oscillator and 
Coulombic spectra of bound states exemplifies the transition from (5) to (7) aX q — 
0. Similarly, we shall not distinguish, at any g > 0, between the wave functions 
pertaining to the symmetric well (5) and to its descendants (7) or (8). 

Returning to the simplest g = models, let us emphasize that they are extremely 
exceptional, possessing 

• all their wave functions in terminating Taylor-series form (note that their fac- 
tors A^^'^\r) are Laguerre polynomials)] 

• all their energies in closed form (note that the HO set forms an equidistant 
family) . 

As a consequence, one should not be surprised by the existence of numerous symme- 
tries (and even supersymmetries [2]) in the underlying Hamiltonians at g = 0. 

At g 7^ 1, many (though not all) of these symmetries become hidden or lost 
(see the monograph [3] for wealth of details). Still, one reveals that the exceptional 
polynomial solutions exist in the form 

N-l 

V'(r) = E h^n^ exp [-BwKB{r)] (9) 

n=0 

at all the finite integers A'^ > 1 and q>l (see the review of this point in our recent 
paper [4]). 
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2 Schrodinger equation at large £ 



For our canonical potential (5), the use of the quasi-exact solution ansatz (9) converts 
the differential equation (1) in algebraic recurrences 
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(10) 



with coefficients 



Cn = (2n + 2) (2n + 2^ + 3), B^ = E - ao (4n + 2^ + 3) 
^ -a^ (4„ + 2e + l) + al- go, A(^^ = -a2 (4n + 2i - 1) + 2aoa^ - gi, 



A^^^ = —aq {An + 2^ + 3 - 2g) + (aottg-i + + ■ ■ ■ + a^-iao) 

n = 0, 1,... . 



(11) 

They form a finite set of algebraic equations which can hardly be solved non- 
numerically at the generic q and A^. In most cases, people only pay attention to 
their very first "square- matrix" special case at g = 1 [3]. 

In what follows, let us admit an arbitrary pair of integers q and and, for 
simplification, accept merely the assumption that the spatial dimension D is very 
large. In the other words, on the basis of the well known formula 



D-3 D-1 D + 1 D + 3 



[12) 



we postulate that these numbers are all very large, i ^ 1. This is a key assumption of 
our forthcoming considerations, inspired by the well known fact that for an?/ potential 
V{r), the practical solution of radial Schrodinger equations is easier in the domain of 
the large angular momenta (a deeper explanation may be found, say, in the randomly 
selected paper [5] or in many other relevant papers with citations listed therein). 
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3 Terminating solutions at large i 



In our present very specific context of tlie incomplete exact solvability, we shouldn't 
be misled by the observation that virtually all the contemporary £ 3> 1 calculations 
are based on the perturbation expansions using the "most natural" artificial expan- 
sion parameter l/£. Rather, we shall follow our older paper [6] (on the q = 2 partial 
solvability) as our most relevant guidance in what follows, having in mind the use of 
a generalized expansion parameter i/fcons* 

In its spirit, our first step will consist in a re-scaling of our over-complete linear 
set (10), Q{E)h = 0, in accord with the simple rule 



2a„ 



(13) 



In this way, all the elements of our non-square band-matrix "Hamiltonian" Q{E) 
become tremendously simphfied in the leading order in D » 1. In effect [4], we then 
have to solve the much easier algebraic problem with N columns and N + q — l rows. 
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where we merely re-scaled the energy E 
linear manner. 
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(14) 



-g-i) and couplings {g^, . . .,gq-2} in 



9k-2 = -otk-iD - 



k-l 



Sfe, 



k 



1,2,. T=(2«+^ 



\ 1/(9+1) 



(15) 



At this stage of development, one does not see any perceivable progress yet. The re- 
quired "multi-spectrum" of q different "multi-eigenvalues Si, . . . , Sg seems obtainable 
only by purely numerical means at all the larger q or N . 
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4 A brief summary of the known non-numerical 
results 



Besides the well known non-numerical q — solutions of the harmonic oscillator, also 
the first nontrivial q = I case need not be discussed too thoroughly. One just solves 
the linear algebraic eigenvalue problem with spectrum which proves equidistant in 
the limit D — > oo (see [7]). For inspiration, let us briefly return to this q = 1 model in 
more detail: In eq. (14), the unknown quantities s = Si represent either the energies 
of the sextic oscillator of eq. (5) or, mutatis mutandis, the charges of the spiked and 
shifted harmonic oscillator [3], the values of which form a finite set. 



We may also very quickly recollect the next q = 2 case where the solution of our 
problem has been found and discussed thoroughly and with direct reference to the 
quartic oscillator potential (7) in 1999 [6]. The first nontrivial form of our equation 
(14) has been solved there in closed form for so many values of N that the results 

could be extrapolated to all = 1, 2, In particular, the resulting energies were 

shown there to form the multiplets 

si = S2 = N - 1, N - 4, N - 7, . . . , -K + 2, N = 2K, 

si^ S2^ N -1,N -A,N -7,...,-K, N^2K + 1, q^2 (17) 

i.e., s ^ N + 2 -3j,j ^1,2,,[(N + 1) /2] at any wave-function degree = 1, 2, . . .. 
In the next step of development, an optimal calculation method has been discovered 
in our subsequent study in 2003 [4]. There, we succeeded in the re- interpretation 
and re-calculation of all the above q = 2 energies as special real roots selected out 
of "hidden-symmetric" complex triplets = s^^^e^^"^^^ , m — 1,2,3. The inter- 
mediate, auxiliary variables s were produced again, numerically, as roots of a set of 
polynomials - 7s^ - 8 = (A^ = 3), s^O - 27s^ + 275^ - 729s = (iV - 4) etc. 

In the same paper, the complete solution of the next problem with g = 3 has been 
offered. The very similar sequence of the secular polynomials F[s'^] has been obtained 
there, with F[s*] = s^-l2s^-Us = at = 3, with F[s*] = si6-68s^2++50625 = 
at A?" = 4 etc. This gives the result 



the presentation of which proves hindered by the occurrence of the other two inde- 
pendent eigenvalues si and S3. Incidentally, the latter quantities coincide and may 
be specified by a closed formula. For our present purposes it is sufficient to elucidate 



si = A^ - 1, iV - 3, iV - 5, . . . - iV + 3, -AT + 1 , q = l. 



(16) 



S2 = Ar-l,A^-5,A^-9,...,A^ + 3- 4 



■A^ + 1" 
. 2 . 



? = 3, 



(18) 
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the A'"— dependence of the resulting multi-spectrum via its first few examples, 
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(20) 



(21) 



etc. As long as wc solved eq. (14) for a sufficiently long scries of "dimensions" N, we 
succeeded in determining the general, extrapolated pattern for s. The paper itself 
should be consulted for more details since the latter spectrum proves to have an 
impressively compact representation in integer arithmetic, with S2 = N + 3 — Aj, 
j = l,2,...,[(iV+l)/2] etc. 

At the time of its derivation, this feature looked ephemeral as definitely failing 
to hold at the next degree q = 4 of the potential. At the same time, although the 
computer-assisted solution of eq. (14) ceased to be feasible, the q = A problem looked 
extremely interesting as long as it involves not only a less appealing polynomial of 
the symmetric well (5) of the eighteen degree, but also much more interesting octic- 
polynomial anharmonic oscillator (7) and, first of all, the phenomenologically most 
important case of the asymptotically cubic force (8). 



5 Brand new result: The case of g = 4 

Before a thorough description of our present continuation of the systematic and 
efficiently computerized symbolic-manipulation study we should re-emphasize that 
at any q >2, our algebraic set of A^ + g — 1 equations (14) is nonlinear. It is formed 
by the sums of the one- and two-term products of the N + q unknown quantities. In 
the latter role we selected the N — 1 arbitrarily normalized Taylor coefficients pj and 
the q multi-eigenvalues Si, S2, . . .Sg. 

In order to convey the feeling of what happens when one chooses the different 
strategics of the usual elimination, we may start experimenting at A^ = 1 an find 
that the only real solution is trivial, Si = S2 = ^3 = S4 = 0. The solution at A^ = 2 
is also unambiguous. Once we abbreviate si = t, S2 = r, S3 = f, S4 = i, we may 
fix the norm by setting pi — 1 and proceed, recurrently, in an upwards direction in 
(14). This gives po = —t while = 1, providing finally the unique real root t — 1 
and, subsequently, full solution with f — r — t — 1. 
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We have seen in ref. [8] that the similar construction is also feasible at TV = 3. 
Proceeding in an upward-downward-symmetric recurent manner we now normalize 
Pi — 1 and infer that po — —1/t while p2 — —1/t. Next we abbreviate tt — ^ and 
re- write the remaining four lines of eq. (14) in the following form, 

rfA = e-2, fi/f^C-^, rt/P^C-^, ft/t^i-2. (22) 

The ratio of the two odd or two even lines eliminates r or f, respectively, and we 
get the same quantity (^ — 2)/(^ — 3). Its next-step elimination gives the desired 
simplification = t^, with the only real solution t = t. Then the first and last line of 
eq. (22) define easily r = r(^) and f = f(^) while, finally, the appropriate insertions 
in one of the middle lines results in the "secular" equation 

p(^t) =t^ -t^ -3t + 2 = (23) 

with the following three real roots, 

h = 2, ^2,3 = ^ (-1 ± v^) (24) 

(cf. also Table 1 below). 

The ehmination of the unknowns becomes almost prohibitively tedious from N — 
5 on. The comparatively high complexity of the (necessarily, computerized) reduction 
of our multi-polynomial problem (14) to the single polynomial "secular" equation 
V{s) = is accompanied by an extremely quick growth of the degree of our secular 
polynomials with N. At the same time, there exists an empirically observed fact [4] 
that, paradoxically, the Grobner-based solution of the next g = 5 problem is in fact 
more easy than its g = 4 predecessor. This underlines the key importance of the 
revealed "missing pattern" in the q — A roots as presented here in Table 1. 

For compensation, the impression produced by the high degree of our polynomials 
V{s) is again strongly weakened when we notice that these functions depend in 
effect just on the powers of the new auxiliary variable z — This has several 

consequences. Firstly, we see that even if all the auxiliary roots z themselves were 
real, the final number of the complex roots s would still be much higher than that of 
their real and, hence, "physically acceptable" partners. Secondly, the formidable task 
of the search for the real roots in the closed form did not prove to be as prohibitively 
difficult as it might have appeared at first sight. 

In Table 1 summarizing the results of our q — A construction, a climax of our 
present effort is perceived in an absolute regularity of all its items. The pattern of 
extrapolation of these results beyond their boundaries set by the computer is already 
fully obvious, 

s = S4(7V) = l (Pjv±V5-g) , 

P(N) = P(iV)o- fe) = 27V + 13 - 5j -10k, Q^j^k) = j - 1 , 

j,k^ 1,2,..., 2j + Ak<N + 5 (25) 
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and does not seem to create any doubts and/or unanswered questions. With respect 
to the non-doubhng of the j — 1 (i.e., Q = 0) roots, the elementary formula 



total # = 



( 



K + 1 
2 
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K 



N + 1 



2 



(26) 



also expresses the total number of the separate items in each column of Table 1, i.e., 
of the real energy roots at each fixed N. 



We reported the progress achieved in the field where the quasi-exact solutions are 
sought for the radial equations where the potentials are "next-to-most-common". 
Our main result is that we were able to construct the energies for the class of the q — A 
models which involves the important and very popular cubic and octic anharmonic 
oscillators. 

Our main task lied in the necessity of making the form of our exact and polynomial 
wave functions ip{r) closed and explicit for a// their integer degrees N = 1, 2, . . .. The 
main difficulty in this direction emerges from the fact that the dimensions of the 
matrices we need (or degrees J\f{N) of the "effective" secular polynomials) seem to 
grow extremely quickly with A^. Unfortunately, we did not find any regularity in the 



series AA(5) = 70, A/'(6) = 126, 7V(7) = 210, J\f{8) = 330, X{9) = 495, ^^{10) = 715 



Our task was quite challenging formally, and we must admit that we did not even 
expect that its solution could appear very soon. Our biggest surprise occurred in the 
form of the explicit factorizability of all the polynomials V{s) over the (sometimes 
called "surdic") field of the quasi-complex numbers a+b with rational coefficients. 

One cannot resist to re-emphasize here that after a certain suitable re-numbering 
and re-grouping of levels, the spectrum of our "solvable" q — A couplings/energies 
remains expressible directly in terms of integers. Such a type of a generalized equidis- 
tance re-emerges also in the next, q = 5 case (i.e., for the class of potentials involving 
the squarc-root-powcr-scrics form of the quartic oscillator, etc). The analysis of g = 5 
already lies beyond the scope of our present study. Even in the purely formal setting, 
it lies on the very boundary of the capacity of the computers and software which 
are at our disposal at present. We were still able to factorizc the corresponding 
effective secular polynomials at a few in ref. [4], and we obtained the regular 
recipe S5 = — 1, iV — 2, — 3, . . . , — A^ + 1 there. One feels how this achievement 
was formidable since at A^ = 7, the extreme coefficient c in the secular polynomial 
F(s) = s^^^ — 600715^^-'^ + . . . + cs possesses as many as 72 decimal digits and, 
hence, looks like a candidate for being placed in the Guiness' book of records in the 
factorization context. 



6 Summary and outlook 



etc. 
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In conclusion, let us point out that our results sample a nice mathematics in 
interplay with a useful physics. Thus, in physics, the equidistance and representation 
of the energies in integer arithmetics in the — > oo limit will enable us to work, in 
any "realistic" dimension D < oc, with perturbation theory without rounding errors. 
In mathematics, the ease of the factorization of polynomials almost certainly reflects 
a hidden symmetry of the Schrodinger equation, but in the light of the nonlinearity 
of its present "algebraization" , we still do not dare to predict any form of its possible 
"explicit manifestation" in the future. 

Besides that "new horizon" , let us also stress once more that our present study has 
been motivated by the disturbing paradox (revealed in [8]) that "phenomenologically 
the simplest" cubic oscillator (such that V{x) ~ a;^ for a; 3> 1) belongs, in terms of 
mathematics, among "the most difficult" examples when its incomplete but exact 
D ^ 1 solvability is concerned. In this sense, we described here a resolution of this 
paradox, showing that the existence of the elementary and exact wave functions ip{x) 
in the large— D regime and for any degree N is admitted not only by the standard 
anharmonic Schrodinger equation with q = 2 (involving the quartic potentials) but 
also by its cubic analogue with q = 4. 
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Table 1. Columns of energy- roots S4 = S4{N) = | (P^ ± \/5 ■ Qj at q — A. 
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